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. Abstract. We present initial results from the M33 Metallicity Project. Out 

' of the thousands of cataloged H II regions in M33, only ~ 30 have electron- 

Cn| ! temperature based abundances in the literature. We have obtained Keck spec- 

q ■ troscopy of a sample of ~ 200 H n regions in M33, with 61 detections of 

' the [O III] A4363 A line that can be used for determining electron tempera- 

, tures, including measurements at small galactocentric radii where auroral lines 

are generally difficult to detect. We find an oxygen abundance gradient of 
—0.027 ± 0.012 dex kpc -1 , in agreement with infrared measurements of the 
neon abundance gradient but much shallower than most previous oxygen gra- 
dient measurements. There is substantial intrinsic scatter of 0.11 dex in the 
(— I ' metallicity at any given radius in M33, which imposes a fundamental limit on 

O |. the accuracy of gradient measurements that rely on small samples of objects. 

I ■ Finally, we present a two-dimensional map of oxygen abundances across the 

southern half of M33 and discuss the evidence for deviations from axisymmetry. 

& ' 

a 

1. Introduction 

>, 

Despite decades of observations of the chemical abundances of stars, H n regions, 
and planetary nebulae, the chemical evolution of galaxies is still a relatively 
data-starved subject. Among disk galaxies (excluding the Milky Way where 
our location makes it difficult to discern the gl obal structure), the best studied 



objec ts are M101 (25 H n region abundanc es; iKe nnicut t. Bresolin. Sz Garnett 
2003 ) and M33 (32 H II r egion abundances; IVilchez et al' 19881 ; Crockett et al 
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enrichment has proceeded. In particular, the mixing of metals through the in- 
terstellar medium (ISM) is not yet very well constrained observationally (e.g., 
Scalo &: Elmegreenl 120041 . and references therein), and whether the spatial dis- 



2006; iMag rini et gill 120071 ) . In both of these galaxies, the existing metallicity 
measurements have a spatial sampling rate of 0.2 kpc -2 . Since H n regions 
and supernova explosions typically affect volumes that are less than a few hun- 
dred pc in radius, these data are too sparsely sampled to reveal how chemical 



tribution of heavy elements is more complex than the typically assumed simple 
exponential gradient with radius is completely unknown. An additional problem 
in the study of extragalactic abundances is the systematic differences between 
the two main techniques for measuring abundances: the so-called "direct" meth- 
ods in which the electron temperature is measured directly using auroral line 
fluxes, and the empirically calibrated strong-line methods that use the flux ra- 
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tios of bright emission lines as a proxy for the oxygen abundance. Motivated 
by these considerations, we have undertaken the M33 Metallicity Project — a 
long term effort to obtain Keck/LRIS spectroscopy of a sample of ;> 500 H n 
regions in M33 that will provide an unprecedentedly large database of uniformly 
determined oxygen, nitrogen, and sulfur abundances in an external galaxy. 



2. Data and Results 



Using the LRIS spectrograph ( Oke et al. 1995) on the Keck I telescope, we have 
observed over 200 H n region s in M33. Details of the obse rvations, data reduc- 
tion, and analysis are given in iRosolowskv &: Simon ( 2007 . hereafter RS07). We 
detected the [O in] A4363 A line and therefore derived direct oxygen abundances 
in 61 of the H n regions, approximately tripling the sample of M33 metallicity 
measurements available in the literature. 

We display our new abundance measurements as a function of radius in 
Figure Q]a. A weak gradient is visible, but the intrinsic scatter in the data 
(which is significantly larger than the observational unce rtainties) dominates 
the pl ot. We therefore must use the method described by lAkritas &: Bershadv 
(1996) to determine the gradient. We measure an oxygen abundance gradient 
of -0.027 ±0. 012 dex kpc" 1 , shallow er than all previous gradient measurements 
except that of Crockett et al. ( 20061 . hereafter C06), and the intrinsic scatter in 
the metallicities at any given radius is 0.11 dex (RS07). 



3. Discussion 



The oxygen abundance gradient that we measure is in good agreement with in- 
frared measurements of the neon abundance gradient bv lWillner &: Nelson-Patel 
(|2002l ) . who found a slope of —0.034 ±0.015 dex kpc -1 . However, our new deter- 
mination does not appear consistent with most previous optical studies of the 
abundance gradient in M33, which generally found steeper slopes. Why does 
our gradient measurement disagree with these other studies? For the objects 
in our sample that have abundance measurements in the literature, our oxygen 
abundances agree within the uncertainties (see RS07), indicating that system- 
atic errors and sample selection biases are probably not to blame. Since we 
primarily observed H n regions in the southern half of the galaxy, and most 
previous studies focused on the northern side, it is possible that a highly asym- 
metric abundance distribution could account for the discrepant results, but this 
hypothesis can only be tested with additional observations. 

We are therefo re left with two possible explanations. The first is that 
IVflchez et al.l ( 19881 . hereafter V88), the largest of the earlier studies, relied on 
photoionization modeling for the abundance of one H n region near the center 
of M33. The high derived abundance for this object significantly steepened the 
overall gradient they derived to —0.10 dex kpc -1 (after rescaling to our assumed 
distance); as V88 noted, the gradient of the outer H n regions in their sample is 
only —0.05 ± 0.01 dex kpc -1 , much closer to more recent measurements. This 
difference suggests that the photoionization model abundances may be system- 
atically higher than direct abundances. Second, the significant intrinsic scatter 
in the abundances revealed by our larger sample prevents measurements based 
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Figure 1. (a) Oxygen abundance gradient in M33. A linear gradient with 
a slope of —0.027 dex kpc -1 is fit to the data (solid line). Regions with sig- 
nificant He ii A4686 A emission are indicated with ope n symbols. The dashed 
and dotted lines represe nt the gradients measured by I Crockett et all {2006) 
and lVflchez et al.1 (|1988l ). respectively. The data are not consistent with the 
steep slopes of ~ —0.1 dex kpc" 1 found by the original M33 studies, (b) 
Monte Carlo simulation showing the distribution of the abundance gradients 
that would be measured for a sample of 10 H II regions drawn randomly from 
our sample of 61. The large width of distribution shows the uncertainty im- 
printed upon the abundance determi natio ns by the underlying variance. The 
measurements of lVfichez et al.l (|1988t ) and lCrockett et all (|2006l ) are indicated 
and are consistent with the distribution given their small sample sizes. 



on small H n region samples from accurately determining the gradient. We plot 
in Figure [U& the outcome of a Monte Carlo simulation in which we measured 
the abundance gradient in M33 by randomly selecting 10 H II regions out of 
our full sample of 61 many times. The broad width of the resulting distribution 
demonstrates that significantly larger samples are needed in order to derive ac- 
curate results. Even with 61 oxygen abundances, we only detect the gradient at 
2.3 a significance! Although our measured gradient is formally o utside the 1 a 
uncer tainties of the gradient slopes reported by V88, C06, and iMagrini et al. 
(120071 ). when the intr insic scatter and small sample sizes are taken into account 
their results are consistent with our data (see Figure [T]6). 

In Figure [21 we display the two-dimensional metallicity distribution we have 
derived for M33 overlaid on an Ha image of the galaxy. As can also be seen in 
Figure [T]a, the highest metallicity H n regions are not located at the center of the 
galaxy, but rather lie along the southern spiral arm at a radius of 1 — 2 kpc. This 
distribution suggests that the material enriched by the most recent generation 
of star formation in the arm has not yet been azimuthally mixed through the 
galaxy. These data represent some of the first evidence for a non-axisymmetric 
abundance distribution in the interstellar medium of M33. As we build up a 
larger data set and cover the northern side of the galaxy, these measurements will 
allow us to constrain the mixing timescale and better understand the processes 
by which heavy elements are transported through the galaxy. 
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Figure 2. Metallicity map of the southwest half of M33. The white circles 
show the locations of our oxygen abundance measurements, and the solid 
black contours represent metallicity. The dotted ellipses are curves of constant 
galactocentric radius (spaced at 1 kpc intervals), and the center of the galaxy 
is marked in the upper left with an "X" . The highest metallicities occur not 
at the center of M33, but in the southern spiral arm. 
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